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In this paper, we studied the semileptonic decays B —)■ D^*h~i2i by using the 
“pQCD+Lattice QCD” method. We made the extrapolation for the six relevant form fac¬ 
tors by using the input values obtained from the pQCD factorization approach in the low 
region of 0 < <7^ < rn^, and the Lattice QCD input at the end-point = g^ax- We 
then calculated the ratios R{D) and R{D*) of the branching ratios B{B —)• and 

found numerically that (1) the “pQCD-i-Lattice QCD” predictions for the branching ratios 
B{B —)■ agree well with the measured values within one standard deviation; and 

(2) the “pQCD-i-Lattice QCD” predictions for the ratios R{D^*'l) are R{D) = O.SSTlgQg® 
and R{D*) = they agree with the data within 2a deviation, in other words, 

one can explain the “R{D^*^ )-puzzle'” in the framework of the standard model. 
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I. INTRODUCTION 

The semileptonic B decays are very important processes in testing the standard model (SM) and 
in searching for the new physics (NP). The B —)■ semileptonic decays have been measured 

by the BaBar in 2012 [1-5], and by the Belle and LHCb Collaborations this year with different 
methods [6-8]. For the ratios of the branching ratios: = B{B -4- D^*^Ti)r) /B{B -4 

D^*Hiyi) with I = (e,/i), the individual measurements [4-8] and the world averages [9] are the 
following: 






0.440 ± 0.072, 

± 0.026, 

0.391 ±0.041 ±0.028, 


BaBar [4, 5], 

Belle [6, 7], 

HFAG average [9], 


0.332 ±0.024 ±0.018, 
0.293;[J;|J^? ± 0.015, 
0.336 ±0.027 ±0.030, 
0.322 ±0.018 ±0.012, 


BaBar [4, 5], 

Belle [6, 7], 

LHCb [8], 

HFAG average [9]. 


( 1 ) 


( 2 ) 


On the theory side, the standard method to study the semileptonic B -4 decays is 

the heavy-quark-effective-theory (HQET). The SM predictions based on the HQET as given in 
Ref. [10] are the following: 


7^(L)) = 0.296 ± 0.016 , n{D*) = 0.252 ± 0.003. 


( 3 ) 
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Two new lattice calculations for the ratio R{D) [11, 12] are of the form 

7^(L)) = 0.299 ±0.011 , 7^(±)) = 0.300 ± 0.008 . (4) 

By comparing the measured values of in Eqs. (1) and (2) with those SM expectations 

in Eqs. (3) and (4), one can see the following points: 

(1) Eor BaBar data only, the measured values of R{D^*^) disagree with the SM predictions in 
Eq. (3) by about 3.4cr [9], this is the so-called “i?(Zl(*))-puzzle”. 

(2) Eor the ratio R{D) and R{D*), the deviation of l.Scr of the Belle data from the HQET 
predictions of R(D^*^) [6, 7] is by itself not significant, but the Belle results do increase 
the tension of the world average values with the HQET predictions slightly. A combination 
for both BaBar and Belle measurements made in Refs. [6, 7] showed that the deviation can 
reach d.Oa. 

(3) The EHCb reported very recently its first measurement for the ratio R{D*) [8], the job 
considered unfeasible at hadron colliders before. The EHCb measurement of R{D*) = 
0.336 ± 0.027 ± 0.030 itself agrees with SM (HQET) at 2.la level, while the combined 
HEAG averages of R{D) and R{D*) as given in Eqs.(l) and (2) do deviate from the SM 
predictions in Eqs. (3) and (4) at 3.9a level [9]. 

How to understand the above -puzzle” ? which has invoked intensive studies in the 

framework of the SM and/or various NP models, for example, in Refs. [13-29]. As is well-known, 
the theoretical predictions for the ratios R{D) and R{D*) have a direct strong dependence on the 
values and the shapes of the B ^ transition form factors F(+,o)(g^), C(g2) and Ao,i, 2 (g^)- 

In Ref. [26], we calculated the ratios R{D^*'>) by employing the pQCD factorization approach, 
where the relevant hadronic form factors are evaluated in the low region by employing the 
pQCD factorization approach and then extrapolated to the higher region by using the pole 
model parametrization [30]. The resulted pQCD predictions are 77(i9) = 0.430^5];o 26 > R{D*) = 
0.301±0.013 [26]. These pQCD predictions agree very well with the BaBar results in Eqs. (1) and 
(2), but we also know that the pQCD predictions for B ^ transition form factors are reliable 
only in the low region [26, 31-34], and the extrapolation from the pQCD predictions at low g^ 
region to the high g^ region by using the pole model parametrization may have large theoretical 
uncertainties. We need something new to improve the reliability of the extrapolation method, the 
Eattice QCD input of F(g^ = g^ax) what we need to determine reliably the size and the shape 
of all the six form factors ±(g^) in the whole range of 0 < g^ < g^ax- 

In Ref. [13], for example, the authors calculated the ratio R{D) by using the hadronic form 
factors computed in lattice QCD in the high g^ region and found R{D) = 0.316 ± 0.012(stat.) ± 
0.007(syst.), which reduced slightly the tension between the SM prediction in Eq. (3) and the 
BaBar result in Eq. (1). 

In this paper, we firstly calculate the form factors F(g^) in the lower g^ < region by using 
the pQCD factorization approach, and then take the F(g^ = g^ax) evaluated in the Eattice QCD as 
additional input to improve the reliability of the pole model extrapolation, and finally to calculate 
the branching ratios B{B —)■ D^*H~i/i) and the ratios R{D^*^). 


II. FORM FACTORS AND THEIR EXTRAPOLATIONS 

In the B meson rest frame, we define the B meson momentum Pi, the momentum P 2 , and 
the polarization vectors e of the D* in the light-cone coordinates in the same way as in Refs. [26, 
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31]. For the B meson wave funetion, we use the same one as being used for example in Ref. [26]: 


= 


\/2iVe 




( 5 ) 


with the distribution amplitude 0 b in the form of 


(pB{x,b) = Nbx‘^{1 


a:)^exp 


■ 

20;^ 


1 

2 


{uJBbf , 


( 6 ) 


where the normalization faetor Nb depends on the values of the parameter ub and deeay eonstant 
/b and is defined through the normalization relation: dx 4 ’b{x, 6 = 0) = /b/(2\/6). We set 

the shape parameter ojb = (0.40 ± 0.04) GeV, in order to estimate the uneertainties of theoretieal 
predietions. 

For the D and D* meson, we use the same wave funetions as being used in Refs. [26, 35] 


$b(p,x) 




V6 


+ mB))0B)(2:), 


\/6 


+mB*)0B*(^) +h{i>D* + mD*)(pD,{x)] , 


( 7 ) 

( 8 ) 


with the distribution amplitudes 


0£){.) (x) 


fpM 

2\/6 


6x{l — x) [1 + Cbi(*) (1 — 2x)] exp 


2 




( 9 ) 


and with the ehoiees of = /b*, 0b* = 0d* = ^d*, and Cn = Cn* = 0.5, uj = 0.1 

GeV as in Refs. [26, 35]. 

In the framework of the SM, the differential deeay widths of the semileptonie B —)■ DWi deeays 
ean be written as [21] 

dr{B ^ DWi) ^ GllKfcp 

dq^ 1927r^m^ \ ) 2q^ 

■{3mf {ml-mlY \Fo{q‘^)\^ 

+ {mf + 2q^)\{q^)\F+{q^)\^y ( 10 ) 


where mi is the mass of the relevant leptons, A(g^) = (m^ + mjj — q^Y — Am\m\) is the phase 
spaee faetor, F(+^o)(?^) are the form faetors. 

For B —)■ D*Wi deeays, the eorresponding differential deeay widths ean be written as [30] 


dFi^{B^ ^ DHvi) 
dq^ 


Gl\V, 


cb\ 


1927i^m% 


1 - 


m 


q. 

2\ a2/„2\ 


" Al/2(g2) 

2g2 

(m| - m|,* - g^)(mB + mB*)^i(g^) 


A(g2 


mB + mr) 


-A2{q^ 


( 11 ) 




















4 


dr±(50 ^ D*Wi) 


GlW A 

1927r3m^ \ 

X |(mf + 2g2) 


2 

Vjq^) 

tub + rriD* 




{rriB + mD*)Ai{q‘^ 

VW) 


( 12 ) 


where A(g^) = (m| + m|,* — g^)^ — is the phase spaee faetor, V (g^) and Aq, 1 , 2 ( 0 '^) are 

the relevant form faetors. The total differential deeay widths ean be defined as 


dr _ d^L dr+ dr_ 

dg2 dg2 dg^ dg^ ’ 


(13) 


The traditional methods to evaluate the form faetors for B D^*'> transitions are the QCD 
sum rule and the Lattiee QCD. Eaeh method has its own advantage region: the QCD sum rule are 
reliable in the low g^ region, but the Lattiee QCD in the high g^ region of g^ ~ 5'max- Refs. [31, 
35], the authors examined the applicability of the pQCD approach to i? —)■ D^*'^ transitions, and 
have shown that the pQCD approach based on the /ct factorization theorem to the B —)■ D^*'^ 
transitions is reliable in the lower g^ region. In Ref. [26, 36], we evaluated the form factors for 
B/Bg ^ ) transitions in the lower g^ region and obtained the pQCD predictions for the 

ratios R{D^*^) and R{Ds*^). 

In order to improve the reliability of the extrapolation of the form factors from the low to high g^ 
region, we here take the Lattice QCD predictions for all six relevant form factors at the end-point 
g^ = ?max the additional input in the fitting process. 

The parametrization of the form factors are rather different in different theories or models. Lor 
B ^ D transition, the form factors _(te) ( or Q(w) ) in the Lattice QCD has simple relations 
with the form factors F+ o(g^) in the QCD sum rule and/or the pQCD approach [14, 21, 27] 


Fo{q^] 


[(1 + r)h+{w) - (1 - r)h_{w)] = Giw), 


2^ 


2^ 


= vr 


1 +w 
1 + r 


h. 


w 


1 


Ue 


h_{w) 


y/mBmp 

rriB - mp 


G{w) A{w), 


(14) 

( 15 ) 


where r = mp/mB, w = {m% + ~ q‘^)/{2mBmp) with g^ = {pb — PnY- The scalar density 

A{w) in Lq. (15) is approximated by a constant value A{w) = 0.46 ± 0.02 [10, 24]. 

Lor the B ^ D* transition, similarly, the form factors L(g^) and 74o,i,2(g^) are related to the 
Lattice QCD ones hA^w) through the following relations [10, 37] 


Ao{q^) = [(1 + w)/iAi(w) - (1 - wr)hA 2 {w) + {r - w)hA^{w )], 


^ 1 ( 7 ^) 

A2{q^) 


1 + r 
1 + r 



(1 +m;)/iai(w), 
{rhA^{w) + hA^{w)), 


( 16 ) 
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where r = rriD^/mB, w = {m\ + m|,* — q'^) / {2mBmB*) with = {pb — Pd*Y, the explieit 
expressions of hv{w) and hA^w) can be found easily in Refs. [10, 37]. 

At the zero-recoil limit: te —)■ 1 and = g^ax’ found simple relations = hMil) = 

/i^ 3 (l) and /iA 2 (l) = 0, and therefore the six form factors at the end-point g^ = q"^^^ , say 
F+,o(gLx)) ^(^Lx) and Ao,i,2(gLx)> depend on two Lattice QCD parameters ^(1) and /iai( 1) 
only. By using the formulae as given in Eqs. (14)—(16) and the following two Lattice QCD 
inputs [14, 37]: 


^(1) = 1.058 ±0.009, /lAi(l) =J^(1) = 0.906 ±0.004 ±0.012, (17) 


one can find the values of the relevant form factors at the end-point g^ = g^ax’ 

^o(gLx) = 0.92 ± 0.02, F+(gLj = 1.21 ± 0.02, V{qiJ = 1.01 ± 0.02, 

.4o(0 = 1.01 ± 0.02, /li(,Lx) = 0.81 ±0,02, 2l2(0 = 101 ± 0,02, (18) 

where ql,„^ = (ms — mu)^ = 11.63 GeV^ for B ^ D transition, q'^^ = (rriB — mD*Y = 10.69 
Gev^ for B ^ D* transition. The uncertainty of the form factors induced by the errors of ^(1) and 
/lAi (1) as given in Eq. (17) is around 0.01 only. We here set conservatively the common error of 
±0.02 for all six form factors in Eq. (18) by taking into account approximately the small variations 
of the central values of ^(1) and /iAi(l) in recent years [11, 14, 37]. 

Now we make the extrapolation from the low g^ to the higher g^ region by using the pole model 
parametrization [30] : 


F{f) 


_ m _ 

1 — a{q^/my) ± b (g^/m|)^ ’ 


( 19 ) 


where F(g^) stands for all the six relevant form factors Fo,+ (g^), C(g^) and Ao,i, 2 (g^). The pa¬ 
rameters a and b in the above equation are determined by the fitting to the values at the seventeen 
points: the first sixteen points in the lower g^ region of g^ < calculated by using the pQCD 
factorization approach, and the last one at the end-point g^ = from the Eattice QCD approach 
in Eq. (18), as illustrated explicitly in Eig. 1. 

In Eig. 1, the dashed curves show the theoretical predictions for the g^-dependence of the six 
relevant form factors in the pQCD approach only [26], while the solid ones denote the six form 
factors evaluated by using the ”pQCD -i- Eattice QCD” method. One can see clearly that: both 
the values and the shapes of all six form factors changed explicitly, in the higher g^ region of 
7 < g^ < gLx’ because of the inclusion of the Eattice QCD input as listed in Eq. (18) at the 
end-point g^ = q^^. 


III. NUMERICAL RESULTS 

In the numerical calculations we use the following input parameters (here masses and decay 
constants in units of GeV) [38]: 

m^o = 1.865, mD+ = 1.870, m^o. = 2.007, mD*+ = 2.010, 

ms = 5.28, mr = 1.777, m^ = 1.35 ± 0.03, fB = 0.21 ± 0.02, fn = 0.223, 

\Vch\ = (39.5 ± 0.8) X 10■^ = 0.287, rs± = 1.638 ps, rso = 1.519 ps, (20) 


while /s* = fD\/mD/mD*, fo* = fos yrncJmD*. 
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(GcV^) B->d‘ (GeV*) B^D* q^ (GeV^) B^D* 

FIG. 1. (Color online) The theoretical predictions for the g^-dependence of the six form factors in the pQCD 
approach (the dashed curves), and the ”pQCD+Lattice QCD” method (solid curves). 


In Fig. 2, furthermore, we show the g^-dependenee of the theoretieal predietions for the differ¬ 
ential deeay rates df/dg^ for B —)■ D*Wi with / = (e“, /i“) and B —)■ D*Ti'r deeays calculated 
by using the pQCD approach only (the dashed curves ) [26], the traditional HQET method [10] 
(the short-dots curve) and the new “pQCD -i- Lattice QCD” method, respectively. 

From the curves in Figs. 1 and 2 one can see that 

(1) For all the six form factors, the g^-dependence of the “pQCD-i-Lattice QCD” predictions are 
similar with those pQCD ones in the lower g^ region, but become rather different in both the 
shape (flat) and the magnitude (smaller) from the previous ones in the large g^ region due 
to the strong suppression effect of the Lattice QCD input at the end-point g^ = g^^^^ in the 
extrapolation. 

(2) For the considered differential decay rates, the “pQCD-i-Lattice QCD” predictions become 
much smaller than the previous pure pQCD predictions in the large g^ region, but approach 
to the HQET ones in the large g^ region. 

In Table I we list the theoretical predictions for the branching ratios of the considered eight 
semileptonic B decays, where the theoretical errors from different sources (say ojb = 0.40 ± 0.04 
GeV, Jb = 0.21 ± 0.02 GeV, |Kfe| = (39.5 ± 0.8) x lO-^ and = 1.35 ± 0.03 GeV) have 
been added in quadrature. As a comparison, we also show the theoretical predictions obtained by 
using the traditional HQET method [10] or the pQCD factorization approach only [26]. The world 
averages as given in PDG 2014 [38] are also listed in last column. One can see from the numerical 
results as shown in Table I that 

(1) For all eight decays, the theoretical prediction from the pQCD approach [26] and the 
“pQCD-i-Lattice QCD” method agree well with each other within one standard deviation. 
They are also well consistent with the world averages as listed in last column [38]. 
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FIG. 2. (Color online) The theoretical predictions for the g^-dependence of dT/dg^ for B —)• D*l~ui 
I = (e, fi) and B D*T~Vr decays in the pQCD approach, the “pQCD + Lattice-QCD” method, and the 
traditional HQET method. Here = 10.69 GeV^ 


TABLE I. The theoretical predictions for the branching ratios (in units of 10 of the B —)• and 

B TVr decays. The previous theoretical predictions in Refs. [10, 26] and the world averages in 

Ref. [38] are also listed as a comparison. 


Channel 


pQCD [26] pQCD+Lattice QCD HQET [10] PDG [38] 



D^T Ur 

n C7+0-34 

o 

00 

1 + 
p p 
to L3 
O CO 

0.57-0.69 

1.03 ±0.22 

B^ 

D+l-ui 

2.03lO;?2 

9 q-j +1.05 
^•'^-'■-0.77 

9 1 Q+0.19 
^•l'i-0.18 

2.19 ±0.12 

B- D^T-Vr 

0.95l°J7 


0.61 - 0.75 

0.77 ±0.25 

B- 

D^l-ui 

9 1 Q+0-99 
"^•-^^-0.76 


2.30 ±0.20 

2.27 ±0.11 

B^ 

D*+T-Ur 

i.seiojs 

-‘-•^^-0.24 

1.21 - 1.35 

1.84 ±0.22 

B^ 

D*+l-ui 

4.52l};^t 

4 70+1-37 

4.94 ±0.15 

4.93 ±0.11 

B- 

D*^T-Ur 

1 47-I-0-43 
-0.40 

1 4n+i-'='7 
1-^11-0.26 

1.31 - 1.48 

1.88 ±0.20 

B- 

D*^l-ui 

-1.41 

c- 1 0+1.49 
^•-*-8-1.28 

5.35 ±0.16 

5.69 ±0.19 


(2) For the branching ratios of the considered semileptonic decays, there is no obvious devi¬ 
ations between the theoretical predictions and the measured ones, since the corresponding 
theoretical errors are still large, around ~ 30 % in size and can not be reduced soon. 


Since the large theoretical errors for the branching ratios are largely canceled in the ratios 
R{D) and R{D*), one can define and calculate these two ratios rather than the branching ratios 
with much better precision. By using the “pQCD -i- Lattice QCD” method we found the following 






















8 


TABLE II. The theoretical predictions for the )-ratios. In last three columns we list the measured 

values and their averages [4, 6, 7]. 

Ratio pQCD[26] pQCD+Lattice HQET [10] HEAG [9] 

n{D) 0.337lg;[]p 0.296 ± 0.016 0.391 ± 0.041 ± 0.028 

n { D *) o.3oil[j:g}^ o.269l[j:g^^ 0.252 ± 0.003 0.322 ± o.ois ± 0.012 



EIG. 3. The 2-dimensionaI R{D) — R{D*) likelihood plane quoted directly from Z.Eigeti’s talk (see foot¬ 
note 1) with the inclusion of the new theoretical predictions obtained by using the pQCD approach (the 
purple box ) or the ”pQCD -i- Eattice QCD” method (the light-blue square) 


theoretical predictions: 


R{D) = 0.337 ± 

_ fl qqy-l-O.OSS 

— U.OO( _o.o37; 

R{D*) = 0.269 ± 0.018(a;B)l|]:[JJ^(/B) 

= 0.269l°;°2i, (21) 

where the remaining major theoretical errors come from the uncertainties of ub = 0.40 ± 0.04 
GeV and fs = 0.21 ± 0.02 GeV. In Table II, we list the theoretical predictions for the R{D^*^)- 
ratios obtained by using the pQCD approach [26], the “pQCD-i-Lattice QCD” method (this work), 
and the traditional HQET method [10]. In last column, we also list the world averaged results as 
given by HFAG [9]. 
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In Fig. 3, we draw the pQCD predietions for both ratios R{D) and R{D*) direetly on the 
2-dimensional R{D) — R(D*) likelihood plane quoted direetly from Z. Ligeti’s talk where 
the purple square and the light-blue square added by us shows the theoretieal predietions with 
the errors obtained by using the pQCD approaeh [26] or the ”pQCD -i- Lattiee QCD” method 
respeetively. 

From Table II and Fig. 3 we find the following points: 

(1) The errors of the pQCD predietions for the ratios R{D^*^) beeome now around 10 %, mueh 

smaller than those for the branehing ratios. But the eurrent 10 % errors are about two times 
larger than the previous one (~ 5 %) as given in Ref. [26], sinee we assigned eonservatively 
an overall error ±0.02 to the Lattiee QCD predietions for the values of the six form faetors 
at the end-point as given in Eq. (18). 

(2) From Fig. 3, one ean see easily that our “pQCD-i-Lattioe QCD” predietions for the ratios 

R{D) and R{D*) do agree with the world averages within 2a level, This faet means that 
one ean explain the ” -puzzle” without the aid from the new physies beyond the SM. 

For the B — D^*Hui deeays, the pQCD predietions do depend on the distribution amplitudes 
of meson as given in Eq. (9). In order to estimate the truneation effeets of the Gegenbauer 
polynomials in Eq. (9), we set the parameter Co = 0 and Co* = 0 for the DAs in Eq.(9), and find 
the new pQCD predictions for the central values of the ratios R{D) and R(D*): 

R{D) = 0.358, R{D*) = 0.282. (22) 

One can see that the pQCD predictions for R{D^*^) become larger slightly than previous ones in 
Eq. (21), ~ 5 % in magnitude and within one standard deviation. It is not hard to understand 
such result: the changes of the form factors and then the branching ratios of the considered decays 
induced by setting = 0 are largely canceled in the ratio R(D^*^). 


IV, SUMMARY 

In this paper, we studied the semileptonic decays B —in the framework of the SM 
by employing the “pQCD-i-Eattice QCD” method. We made the extrapolation for the six relevant 
form factors ^(+,0) (g^), C(g^) and Aq, 1 , 2 ( 0 '^) by using the input values at sixteen points from the 
pQCD factorization approach in the lower region ( 0 < < m^), and the ones at the end-point 

g^ = g^a^ from the Eattice QCD predictions as given in Eq. (18). We then calculated the branching 
ratios B{B —)■ and the ratio R{D) and R{D*) by employing the new ”pQCD-i-Eattice 

QCD” method. 

Erom the numerical results and phenomenological analysis we found the following points: 

(1) By using the Eattice QCD predictions for the six form factors as the input at the end-point 
g^ = g^ax’ *^be reliability of the extrapolation of the form factors from the low to the high g^ 
region are improved effectively. 

(2) The “pQCD-i-Eattice QCD” predictions for the branching ratios B{B —agree 
well with the measured values within one standard deviation. 


* This is the figure appeared in Z. Ligeti’s talk given at the ICTP summer school on particle physics, June 15-26, 
Trieste, Italy, and where the HQET predictions for the relevant experimental measurements from BaBar, 

Belle and LHCb collaborations, and the world averages [9] have been illustrated explicitly. 
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(3) The ”pQCD+Lattice QCD” predictions for the ratios are the following 

R{D) = 0.33710:°38, R{D*) = 0.269”. (23) 

They agree with the world averages of those currently available experimental measurements 
within 2cr deviation, in other words, one can explain the so-called ”R{D^*^)-puzzle” in the 
framework of the standard model. 
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